Combinatorial Search for Optimal Hydrogen-Storage Nanomaterials Based on 

Polymers 
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We perform an extensive combinatorial search for optimal nanostructured hydrogen storage mate- 
rials among various metal-decorated polymers using first-principles density-functional calculations. 
We take into account the zero-point vibration as well as the pressure- and temperature-dependent 
adsorption-desorption probability of hydrogen molecules. An optimal material we identify is Ti- 
decorated czs-polyacetylene with reversibly usable gravimetric and volumetric density of 7.6 weight 
percent and 63 kg/m 3 respectively near ambient conditions. We also propose "thermodynamically 
usable hydrogen capacity" as a criterion for comparing different storage materials. 

PACS numbers: 68.43.Bc, 71.15.Nc 



Hydrogen storage is a crucial technology to the devel- 
opment of the hydrogen fuel-cell powered vehicles Q, @ • 
Recently, nanostructured materials receive special at- 
tention because of potentially large storage capacity 
(high gravimetric and volumetric density), safety (solid- 
state storage), and fast filling and delivering from the 
fuel tank (short molecular adsorption and desorption 
time) HS@. 

However, when the thermodynamic be- 
havior of the gas under realistic environments is taken 
into account, the usable amount of hydrogen with these 
nanomaterials falls far short of the desired capacity for 
practical applications and search for novel storage ma- 
terials continues worldwide 0, S @ • It is to be em- 
phasized that hydrogen storage in nanostructured mate- 
rials utilizes the adsorption of hydrogen molecules on the 
host materials and its thermodynamic analysis is distinct 
from that of metal or chemical hydrides. Each adsorption 
site on the nanomaterial behaves more or less indepen- 
dently and the probability of the hydrogen adsorption 
follows the equilibrium statistics which is a smooth func- 
tion of the pressure and temperature. There is no sharp 
thermodynamic phase transition between the gas and the 
adsorbed state of H2 , in contrast to the case of metal or 
chemical hydrides where an abrupt phase transition oc- 
curs at well-defined pressure at a given temperature 01 • 
With this caveat, a general formalism applicable to the 
hydrogen adsorption on nanomaterials was derived in the 
present study from the grand partition function with the 
chemical potential determined by that of the surround- 
ing H2 gas acting as a thermal reservoir. As each site 
can adsorb more than one H2 molecule, information on 
the multiple adsorption energy is necessary. (The situ- 
ation is analogous to the O2 adsorption and desorption 
on hemoglobin which can bind up to 4 O2 molecules.) 
In equilibrium of the H 2 molecules between the adsorbed 
and desorbed (gas) states, the occupation (adsorption) 
number f is obtained from f—kTdlnZ/dfi, where Z is the 
grand partition function, fj, is the chemical potential of 
H2 in the gas phase at given pressure p and temperature 



T, and k is the Boltzmann constant. Here, f per site is 
reduced to 



f 



(1) 



where £1 is the adsorption energy per H2 molecule when 
the number of adsorbed molecules is / and gi is the multi- 
plicity (degeneracy) of the configuration for given I. The 
summation is over all different configurations up to the 
maximum number (N max ) of adsorbed molecules. 

Another important thermodynamic feature in the H2 
adsorption energetics is the zero-point vibrations of the 
H2 molecules with respect to the host metal atom (e.g. 
Ti) on which H2's sit. (The zero-point vibration within 
the H2 molecule, on the other hand, exists in both gas 
and adsorbed states and cancels out in the calculation 
of f.) The actual energy £; to be used in Eq. (1) is the 
static adsorption energy (usually calculated and reported 
in the literature) minus the zero-point vibration energy 
which sums up to as large as 25 % of the static adsorption 
energy according to our calculation, a value not to be 
neglected at all. 

Considering these thermodynamic aspects, we paid at- 
tention to the fact that simple polymers decorated with 
light transition metal atoms may be superior to other re- 
cently reported nanomaterials such as Ti-decorated nan- 
otubes or Sc-decorated fullerenes 0] in terms of 
usable gravimetric and volumetric density. The basi- 
cally one-dimensional nature of polymers is advantageous 
for compact storage, with a very small number of car- 
bon atoms needed to accommodate a decorating metal 
atom which attracts hydrogen molecules. Furthermore, 
entangled chains of long polymers form a solid struc- 
ture ideal for safe handling of the hydrogen. A sys- 
tematic approach was employed to search for optimized 
high-capacity hydrogen storage nanostructures based on 
polymers. For the supporting backbone materials, we 
first considered trans- and cis-polyacetylene (among lin- 
ear carbon chains), polyaniline, polyphenol, poly para 



phenylene, and poly ether ether ketone (chains of hexag- 
onal rings) , and polypyrrole and polythiophene (chains of 
pentagonal rings). For decorating transition metals, we 
initially chose all light transition metal elements starting 
from Sc in the periodic table. Various possible adsorption 
sites of the transition metal atoms were tested for each 
case. The maximum number of adsorbed H2 molecules 
also varied (up to six) at different sites. In short, the 
total combinatorial number in our study exceeded one 
thousand. In practice, we were able to reduce the number 
considerably by eliminating obviously unfavorable cases 
using a few test calculations of the adsorption energy 
and structural stability. Many kinds of pentagonal and 
hexagonal ring chains were ruled out. For decorating 
atoms, only Sc, Ti, and V atoms passed the first-round 
candidate screening test. Such a combinatorial search 
for the optimized material and geometry yielded a few 
promising nanostructures for hydrogen storage. 

We employed spin-polarized first-principles electronic 
structure calculations based on the density-functional 
theory |l3|. The plane-wave based total energy mini- 
mization |l4| with the Vanderbilt ultrasoft pseudopoten- 
tial 0| was performed. The generalized gradient ap- 
proximation (GGA) of Perdew, Burke, and Ernzer- 
hof (PBE) [13 was used in the calculations. The kinetic 
energy and the relaxation force cutoff were 35 Ry and 
0.001 Ry/a.u., respectively. For periodic supercell calcu- 
lations, the distance between polymers was maintained 
over 10 A in all cases. 

The best candidate material we found in our search 
using the total energy calculations was czs-polyacetylene 
decorated with Ti atoms whose structure after the H 2 
molecule adsorption is presented in Fig. 1(a). The 
binding energy of a Ti atom on this polymer is 2.4 eV. 
The structure has about 2 wt% higher storage capac- 
ity than Ti-decorated trans- polyacetylene. (Since as- 
synthesized polyacetylene is of cis-type, it is in princi- 
ple possible to attach Ti atoms to cis-polyacetylene al- 
though trans- polyacetylene is a more stable structure.) 
N max for this structure is five and the H 2 molecules 
are compactly adsorbed on both sides of the polyacety- 
lene plane. The molecular formula corresponding to this 
structure is (C4H4-2Ti-10H2) n . The maximum gravi- 
metric density (G max ) is defined by the weight ratio of 
10H 2 to C 4 H 4 -2TM0H 2 , which is 12 wt% as shown in 
Table I. G max for other materials is calculated in the 
same way. We also present other important polymer ge- 
ometries in Fig. 1 with the maximum number of H2 
molecules attached to the decorating Ti atoms. The cal- 
culated (static) adsorption energies per H2 as a func- 
tion of the adsorption number are presented in Fig. 2 
for easy comparison among different materials. In cis- 
polyacetylene, for example, they are 0.55, 0.58, 0.48, 0.42, 
and 0.46 eV/H 2 for I =1, 2, 3, 4, and 5, respectively. As 
pointed out in previous works, the adsorption of a large 
number of H2 molecules presumably occurs through the 
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FIG. 1: (color online) Atomic structures of the Ti-decorated 
polymers with the maximum number of H2 molecules at- 
tached to Ti atoms. Green, blue, purple, yellow, and red 
dots indicate the carbon atom, titanium atom, nitrogen atom, 
hydrogen atom composing the polymer, and the molecular 
hydrogen, respectively, (a) cis-polyacetylene with five H2 
molecules attached per Ti atom. Hys are shown on both 
sides of the (somewhat distorted) polyacetylene plane. In the 
rest (b)-(f), Hys are shown only on one side of the polymer for 
visual clarity, (b) irans-polyacetylene with Ti atoms located 
out of the plane of the polymer chain, (c) trans-polyacetylene 
with Ti atoms in the plane of the chain, (d) polypyrrole with 
Ti atoms out of the pentagonal plane, (e) polyaniline with Ti 
atoms in the hexagonal plane, (f) polyaniline with Ti atoms 
out of the hexagonal plane. 

Dewar-Chatt-Duncanson coordination or Kubas interac- 
tion (lilll^.l2fj| . We found the elongation of H2 molecules 
by ~T0 % through electron back donation from metal d 
orbitals to the antibonding hydrogen s orbitals, which 
supports these theories. 

We chose to present in Fig. 2 the static adsorption en- 
ergy following the usual practice in the literature 0, . 
After subtracting zero-point vibration energies (25 % of 
the static adsorption energy) for all structures, we ob- 
tained the true dynamic adsorption energy (ei) to be 
used in Eq. (1). For instance, the zero-point vibration 
energy per H2 molecule for cis-polyacetylene was 0.09 
eV for H2 on top of the Ti atom and 0.12 eV for H2 at- 
tached to the side [2l| . We employed the experimental 
chemical potential in the literature |22[ in the calcula- 
tion of f. The degeneracy factor was approximated by 
the number of calculated local energy minima for given 
I. The largest gi we found was 3 and, since the expo- 
nential factor e'( M ~ £i )/ fcT dominated, gi's turned out to 
give a minor correction to the result. The occupation 
number f as a function of p and T for representative 
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FIG. 2: (color online) Calculated static adsorption (bind- 
ing) energy per H2 molecule for polymers decorated with Sc, 
Ti, or V atoms. The average binding energy per H2 is plot- 
ted up to the maximum number of adsorbed HVs allowed for 
each species. PA, PPY, and PANI stand for polyacetylene, 
polypyrrole, and polyaniline, respectively, -out and -in mean 
out-of-plane and in-plane configurations as previously shown 
in Fig. 1, respectively. 



nanomaterials is presented in Fig. 3. The occupation- 
pressure-temperature (f-p- T) diagram of the nanomate- 
rial storage in Fig. 3 is the counterpart of the widely- used 
pressure-composition isotherms (PCI) in metal hydride 
storage To obtain the usable amount of hydrogen, 
it is necessary to specify p and T at the time of adsorp- 
tion (filling) and desorption (delivering from the storage 
tank) . Since an internationally agreed-upon standard has 
not been set up, we propose to use the adsorption con- 
dition of 30 atm and 25 °C and the desorption condition 
of 2 atm and 100 °C, abbreviated to 30-25/2-100. These 
numbers, which may be revised in the future by consen- 
sus, are based on information in the literature |U 0] 
and reflect practical situations in gas filling and vehi- 
cles operations. 30 atm for adsorption and 1.5 atm for 
desorption were used in Ref. 24, but they did not take 
advantage of the temperature variation. 

We adopted an easily achievable temperature range of 
25-100 °C here. Then, f at the condition of 30-25 minus 
f at 2-100 is the available number of H2 molecules per 
site. These numbers are listed in Table 1. For compar- 
ison, the same numbers for the Sc-decorated fullerene 
(C4SB 12SC 12) [l2| and the Ti-decorated carbon nan- 
otube [11| are presented as well. We confirm that our 
results for the static adsorption energy of these materials 
agree with reported values 0, 0] . When converted to 
the gravimetric density, Ti-decorated cis-polyacetylene 
stores usable H2 molecules of 7.6 wt% out of the maxi- 
mum density of 12 wt%, which is much greater than, say, 
the goal of 6 wt% by the year of 2010 set by the Depart- 
ment of Energy (DOE) of US H^. The Ti-decorated cis- 
polyacetylene is the best candidate material for hydrogen 



TABLE I: Hydrogen storage capacity of representative nano- 
materials from GGA calculations. PA, polyacetylene; PPY, 
polypyrrole; PANI, polyaniline; CNT, carbon nanotube. All 
are decorated with Ti except for Sc-decorated C48B12. -out 
means an out-of-plane configuration described in Fig. 1. N a d s 
and Ndes are the numbers of attached HVs per site at the con- 
dition of adsorption (30 atm-25 °C) and desorption (2 atm- 
100 °C), respectively. N usc is the practically usable number 
(N a( j s — Ndes) and N max is the maximum number of adsorbed 
HVs. G and V are gravimetric and volumetric density, re- 
spectively. 



Materials 






Guse/Gmax 


Vuse /Vmax 








(wt%) 


(kg/m 3 ) 


cis-PA 


5.00-1.84 


3.16/5 


7.6/12 


63/100 


PPY 


3.00-0.05 


2.95/3 


4.9/5 


33/34 


PANI-out 


3.00-0.96 


2.04/3 


4.1/6 


31/46 


C48B12SC12 


2.68-0.02 


2.66/5 


4.7/8.8 


23/43 


CNT 


1.95-0.35 


1.60/3 


4.1/7.7 


not available 



storage in the l ist, better than that of Ti-decorated car- 
bon nanotubes jllj or Sc-decorated fullerenes [l2| ■ Note 
that 60 % desorption of H2 is achieved here at a tem- 
perature as low as 100 °C, which is considerably lower 
than the dissociation temperature of usual metal hy- 
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FIG. 3: Occupation number-pressure-temperature (f-p-T) di- 
agram of the hydrogen storage in representative nanomate- 
rials. The ranges of the pressure and the temperature cover 
typical conditions of filling and delivering from the storage 
tank, (a) Ti-decorated cis-polyacetylene. (b) Ti-decorated 
polypyrrole. (c) Sc-decorated C48B12. (d) Ti-decorated car- 
bon nanotube. 
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drides. If we were to slightly raise the desorption tem- 
perature to 130 °C, the usable gravimetric capacity of 
cis-polyacetylene would reach 9 wt%. 

The volumetric density of the hydrogen storage is dif- 
ficult to evaluate and not usually reported in the litera- 
ture. To estimate the volumetric density, we assume that 
a hydrogen molecule adsorbed at the top of one unit cell 
and another adsorbed at the bottom of the next unit cell 
are separated by a van der Waals distance (~3.4 A). The 
calculated usable volumetric density is 63 kg/m 3 , which 
is higher than the 2010 goal of 45 kg/m 3 set by the DOE 
of US 0. 

So far, we have demonstrated an enormous potential 
of polymers as a hydrogen storage medium. A serious 
problem yet to be overcome in practice is the attack of 
oxygen or other ambient gases. Polyacetylene is known to 
be especially vulnerable to it. One possible morphology 
to avoid oxidation is a dense matrix of polymer which 
is made permeable to H2, but not to O2. Clustering 
of decorating metal atoms is another obstacle in mate- 
rial fabrication [25| . In our simulation which allowed for 
the relaxation of atomic positions, individually dispersed 
metal atoms did show local stability in linear polyani- 
line. However, how to suppress the aggregation of metal 
atoms in polymer matrices in general is a difficult ques- 
tion to be answered. As a brief summary of the experi- 
mental situation, we want to point out that the hydrogen 
storage in polymer-dispersed metal hydrides was studied 
before pif. It was hoped that the storage capacity of 
metal hydrides might be enhanced by incorporating a 
low-density polymer that could interact with the hydride 
on a molecular level and store additional hydrogen within 
the polymer structure. Although increase in the hydro- 
gen release was found, the maximum amount reached 
only 0.36 wt% in experiment. We believe that the low 
capacity is probably due to the poor morphology, i.e., too 
large (>10 fj,m) Ti-polymer particles. 

In summary, we have carried out first-principles elec- 
tronic structure calculations for hydrogen binding to 
metal-decorated polymers of many different kinds. When 
the thermodynamic behavior of hydrogen molecules un- 
der realistic conditions is considered, the Ti-decorated 
cis-polyacetylene is found to have the highest usable 
gravimetric and volumetric density among nanostruc- 
tures reported so far. We also propose the l-p-T diagram 
as a criterion for evaluating usable capacity at ambient 
conditions. It remains to be a challenge for experimental- 
ists to fabricate a structure of individually dispersed Ti 
atoms on polymer as much as possible in order to achieve 
significantly improved storage capacity. 
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of MOST/KOSEF and the Korea Research Foundation 
Grant No. KRF-2005-070-C00041. Computations are 
performed through the support of KISTI. 
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